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Abstract: This survey treats the problem of ruin in a risk model when 
assets earn investment income. In addition to a general presentation of 
the problem, topics covered are a presentation of the relevant integro- 
differential equations, exact and numerical solutions, asymptotic results, 
bounds on the ruin probability and also the possibility of minimizing the 
ruin probability by investment and possibly reinsurance control. The main 
emphasis is on continuous time models, but discrete time models are also 
covered. A fairly extensive list of references is provided, particularly of pa- 
pers published after 1998. For more references to papers published before 
that, the reader can consult [47]. 
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1. Introduction 



The problem of ruin has a long history in risk theory, going back to Lundberg 
[40]. In Lundberg's model, the company did not earn any investment on its 
capital. An obvious reason for this assumption, although there may be other 
reasons as well, is that the mathematics is easier, and back in the first half of 

"This survey was written while the author was a visiting professor at the departments of 
mathematics and statistics at The University of Chicago 
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the 20th century the theory of stochastic processes was far less developed, and 
also far less known, than it is today. The first attempt to incorporate investment 
incomes was undertaken by Segerdahl in [02]. Segerdahls assumption was that 
capital earns interest at a fixed rate r. This model was further elaborated in 
[22, 28], and in a somewhat more general form in [18]. The books [4, 54] both 
have sections devoted to it, and it remains very popular even today. We will 
repeatedly return to it in this survey. 

Inspired by ideas from mathematical finance, in [45] a model was suggested 
where capital is allowed to be invested in risky assets. Our starting point here 
will be a slightly restricted version of this model. 

For a survey of the theory before 1998, the reader is referred to [47]. Since 
1998, there are particularly three new developments that has influenced and 
given new vitality to this topic. 

1. The emphasis on heavy tailed claim distributions. 

2. The Gerber-Shiu penalty function. 

3. The possibility to influence the ruin probability by control of the risky 
investments and possibly reinsurance. 

In 1998, the only papers dealing with these items in the context of this survey 
were [3, 7, 35]. 

In order not to be bogged down in technicalities, the reader is referred to the 
references for special cases and detailed assumptions. 

To make the ideas transparent, we introduce the risk process by means of 
two basic processes, i.e. 

• A basic risk process P with Pq = 0. 

• A return on investment generating process R with Rq = 0. 

It is assumed throughout in this survey that P and R are independent. If P and 
R belong to the rather general class of semimartingalcs, then under some weak 
additional assumptions we can define the risk process as 



where R = \og£{R) is the logarithm of the Doleans-Dade exponential of R, 
i.e. Vt = £(R)t, satisfies the stochastic differential equation dV t = Vt-dRt and 
Vo = 1. 

In this survey we shall typically assume that P and R are of the forms 




(1.1) 



so that Yq = y. By [47] the solution of this equation is given as 




Pt 




(1.2) 



i=i 



Ri 



rt + a R W Rit , 



(1.3) 
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where Wp and Wr are Brownian motions, N a Poisson process with rate A and 
the {Si} are positive, independent and identically distributed random variables 
(i.i.d.) with distribution function F. Furthermore, Wp, Wp, N the {Si} are all 
independent. The idea is that p is the premium rate, the {Si} are claims while 
Wp represents fluctuations in premium income and maybe also small claims. 
The return on investment generating process R is the standard Black Scholcs 
return process. With these assumptions, Y becomes a homogeneous, strong 
Markov process, a fact that allows us to draw on the vast literature on Markov 
processes. When R follows (1.3), R\ = Rt — = {r — 5 cj^i + apWn.t- 

A few papers [45, 49, 75, 76] generalize the return on investment process R 
to the jump-diffusion 

N R ,t 

Rt = rt + a R W R , t +J2 S ^ ( L4 ) 
i=i 

where X)i=a* &R,i 1S another independent compound Poisson process. Letting 
Fp(x) = P(Sp < x), in order to avoid certain ruin caused by losing everything 
in an investment shock, it is assumed that Fp(—1) = 0. At an even higher level 
of generality, in [46, 48], P and R are independent Levy processes. 

The time of ruin is defined as T = mi{t : Y t < 0}, with T — 00 if Y stays 
positive. The probability of ruin in finite versus infinite time is then defined as 

i){t,y)=P{T<t\Y =y) and ij)(y) = P(T < oo\Y = y). 

Mathematically, ip{y) is the easiest, and probably as a consequence it has by far 
been the most popular in the literature. 

A more general concept that also allows for problems relating to the time of 
ruin, the size of the deficit at ruin or the surplus immediately before ruin, is the 
Gerber-Shiu penalty function [2'4]. It is given as 

= E[g{Y T _, \Y T \) e - aT l {T<oo} \Y = y], 

where g is a nonnegative function and a > 0. Various choices of g and a allow 
for the computation of many interesting quantities related to ruin [9] . Note that 
with .g(yi,y 2 ) = 1 and a = 0, $o(y) = ip(y)- Another choice is 5(2/1,2/2) = 
l{ vi <u,y 2 <v} and a = giving 

H(y, u, v) = P(T < 00, F T - < u, \Y T \ < v\Y = y). (1.5) 

Papers with material on $ a {y) or special cases of it are [9, 12, 70, 71, 72, 73, 
74, 75, 76, 77]. 

A somewhat different ruin concept introduced in [19] and [57] is the model 
for absolute ruin. Let y > and 

Y t = y + P t + f Y„-l {Y ,_ <0} dRi, a + [ (F s _ - //!l,v .„:<//>'.,.,. (1-6) 
Jo Jo 
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where Ri s t = r±t and i?3 jt = r-^t + upWp.t. The idea is that when capital is 
negative, money can be borrowed at rate T\. When capital is between and y, 
it is kept in the company without earning any investment income, while excess 
capital over y are invested in a possibly risky market. If P follows (1.2) with 
op = and for some T A , Y t a < —p/ri, premium income is no larger than 
interest on debt, and consequently Y drifts towards minus infinity. The time 
T A is called the time of absolute ruin, and ?p A (y) = P(T A < 00) the probability 
of absolute ruin. When op > the concept of absolute ruin is less meaningful 
since whatever small Y t is, it can with some positive probability drift back into 
positive values again. 

2. Some general results in the Markov model 

Unless otherwise stated, it is assumed that P and R follow (1.2) and (1.3). 
Then it follows from [46, 52] that under weak assumptions, ip(y) < 1 if r > \o R 
(and in fact when cr R > then ip(y) < 1 if and only if r > \cr R ). In this case, 
under weak assumptions, ip is twice continuously diffcrcntiable on (0, 00) and is 
a solution of the equation, see [25, 49, 69] and in particular [26], 

my) = -XF(y), (2.1) 

with boundary conditions 

lim i/j(y) = and tp(0) = 1 if a P > 0. 

y^oo 

Here L is the integro-differcntial operator 

Lh(v) = \{°p + °W)h"(y) + (P + ry)h'(y) + A j* h(y - x)dF(x) - Xh(y), 

and F(y) = 1 — F(y). Sometimes it is more convenient to work with the survival 
probability 4>{y) = 1 — ^{y), in which case (2.1) becomes 

I4(y) = 0. (2.2) 

When R instead follows (1.4), it was shown in [75] that under rather strong as- 
sumptions, "0 is twice continuously diffcrcntiable and satisfies Lpip(y) = —XF(y), 
where 

/oo 
^ h(y(l + x))dF R (x)-X R h(y). 

Equations similar to (2.1) also hold for other relevant quantities. An example 
is the decomposition of the ruin probability into ip{y) = ^diy) + ^s{y), where 
4>d{y) is the probability that Y will become negative due to a drift in Wp, while 
ip s {y) is the same, but due to a claim, i.e. a jump. Under suitable differentiability 
assumptions, it is shown in [11] that ipd satisfies Lipdiy) = with the same 
boundary conditions as ip, while ip s satisfies (2.1), but with ip s (Q) = 0. 
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A further example is the calculation of the Gerber-Shiu penalty function. 
With up = and some rather strong assumptions on the distribution F and 
the function g, it is shown in [9] that $> a (y) satisfies 

/•OO 

L* a (y) - = -A / g(y,x - y)dF(x), (2.3) 

Jy 

with boundary condition lim^oo $ Q (y) = 0. The extension to ap > seems 
rather straightforward, in which case the boundary conditon <1> Q (0) = <?(0, 0) 
must be added. Also, using methods such as in [26] , the assumptions can prob- 
ably be weakened considerably 

In the absolute ruin problem there are three equations Li^pf-(y) = — \F(y), 
where the operators Li are as L above, but referring to the different Ri where 
R-2.t = 0. The solution can then be found by imposing proper boundary and 
continuity conditions. 

Following ideas from [28], it was shown in [45] that the ruin probability can 
be written as 

*<») = E w H mT«»Y {2A) 

where H is the distribution function of the perpetuity 
X= / e - {r -^ t+,jRW ^dP t . 



Actually, (2.4) holds whenever P and R are independent Levy processes, see [46]. 

Corresponding to (2.1), the finite time ruin probability ip(t,y) should be the 
solution of the partial integro-diffcrcntial equation 

- Q- t Ht, y) + L4,{t, y) = -XF(y), (2.5) 

with the additional boundary condition ip(0,y) = l{ y<0 }. Here the operator L 
acts on the y variable. When ap = ap : = standard methods can be used to 
prove that the ruin probability is sufficiently smooth and satisfies (2.5), but as for 
the infinite time case this problem is much more complicated when ap + a R > 0, 
and for this case there are no results as of the existence of a smooth solution of 
(2.5). In [42] an equation similar to (2.5) for diffusions is discussed. 

3. A discrete time model 

With {tj} the jump-times of P and To = 0, setting X n = Y Tn yields 

X n = eA»-«^-i |x„_x + jf " e-^- A ^-,) dP \ = £ n x„_! + Tj n . (3.1) 

Here the sequence {(£ n , rjn)} is i.i.d., but ^ n and r\ n are themselves not indepen- 
dent. However, {^niVn) is independent of X n -\. 
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Assume that ap = 0, in which case 

ip(y) = P(X n < for some n), 

hence we can alternatively focus on the sequence {A„}. An advantage with this 
setup is that the claim-number process N does not have to be a Poisson process, 
any renewal process will do, [8, 77]. Actually, (3.1) can serve as a link between 
the time continuous and the time discrete frameworks. To see how, consider the 
time discrete process {X n } where X n is total capital at time n. Let {(R n ,P n )} 
be i.i.d. where R n is the interest rate between time n — 1 and n, and P n is 
premium minus claims in the same time interval. Then 

X n = (1 + FQXn-t + P n = C„X„_! + T) n (3.2) 

with £„ = 1 + R n and r\ n = P n . A major difference between the time-continuous 
and the time-discrete setup is that in the time-discrete case it is common to 
assume that R n and P n are independent, i.e. that and r\ n are independent. 
The rationale for this assumption is that premiums earned and claims paid do 
not influence investment income until the beginning of the next period. This 
is the assumption in [43, 64, 65], but not in [44] where a more general setup 
was used in order to deal with the continuous time case as well, including when 
ap 7^ 0. Also, in [16], the independence assumption of the is relaxed. In [74] 
a discrete time version of (2.4) is given, while [16, 43, 44, 64, 65] all focus on 
the asymptotics of ruin probabilities, to be briefly reviewed below. 

4. Analytical and numerical solutions 

Since most analytical results are for the infinite time horizon, we start with this 
case. When A = the process is a pure diffusion process and the problem is 
solvable. A solution can be found in [45] or [55]. When A > analytical solutions 
can only be obtained in a few rather simple cases, and most of these solutions 
are very complex. In all known solutions it is assumed that ap = 0, so in the 
sequel we shall therefore tacitly let ap = 0. 

We have already mentioned Scgerdahl's classical work when ap = and 
exponentially distributed claims. In [49] these solutions were extended to the 
case when claims are mixtures of two exponential distributions, as well as to the 
case when they are Erlang(2) distributed, i.e. they can be represented as a sum 
of two independent exponentials. Extensions beyond that seem very difficult 
though. In [49] the case with ap > and claims exponentially distributed was 
also solved, and this solution was extended in [10] with separate solutions for 
i/j d (y) and tp s (y). 

Again with ap = 0, the absolute ruin problem when ap = and exponential 
claims was solved in [19]. This result was extended to ap > and y = in [24]. 
Another extension can be found in [12] where rather explicit expressions for the 
Gerbcr-Shiu penalty function are given for the case with y — oo and ap = 
and claims exponentially distributed. Similar results for the case (1.5) are found 
in [73]. 
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In the finite time horizon case analytical solutions are hard to come by. Using 
(2.5), in [36] the Laplace transform of the survival probability cj>(t, y) = l—ip(t, y) 
is found when ap = ap = and claims are exponentially distributed. However, 
this transform involves a ratio of confluent hypergeometric functions, and is 
therefore difficult to invert, the exception is when A = r in which case the 
solution is rather simple. A different method for the same problem is used in 
[1] who provide a recursion for <p(t, y) when A = kr for some positive integer k. 
This recursion is solved and exact solutions given for k = 1 and k = 2. 

The value of continuous time ruin theory is mostly due to its simplicity as a 
concept together with its complexity as a mathematical problem, and this can 
explain why the issue of computing numerical values has received comparatively 
little attention. In [51], following an idea from [63], but allowing a P > 0, using 
integration by parts the equation (2.2) was turned into a Volterra integral equa- 
tion and methods from numerical analysis was used to solve this numerically. 
In the finite time case several methods have been proposed when ap = ap = 0, 
see e.g. [6, 13], and for a somewhat more general discrete time model [17]. These 
methods are rather intuitive in nature and not based on any particular known 
procedure from numerical analysis. Their efficiency arc therefore somewhat low, 
but as a bonus they provide upper and lower bounds. 

An alternative to traditional numerical methods is Monte Carlo simulation, 
which is particularly well suited for finite time ruin problems. For infinite time 
ruin problems some care has to be taken as to when the simulation should stop. 
An alternative is to simulate under an equivalent measure P so that P(T < 

oo\Yq = y) = 1. Then with M t equal to restricted to a{Y s : s < i}, 

ij(y)=E[M^\Y =y], 

hence what is required is repetitive simulation of Mt under P. The catch is that 
this kind of importance sampling is very dangerous, and can lead to completely 
wrong results, see e.g. [50] for examples. In [5] a somewhat complicated but 
efficient change of measure is given for the case with ap — ap = and light- 
tailed claims. The change of measure method also works for finite time ruin 
probabilities, in which case ip{t,y) = E[M^ At \Yo = y]. 

5. Asymptotic results 

If there is only a few papers dealing with numerical solutions, there are certainly 
a fair amount dealing with asymptotic results, i.e. the behaviour of ip(y) or 
ip(t, y) when y gets large. To present a few of the most prominent or recent 
results we shall need some definitions. As before let F{x) = 1 — F(x) and let 
F* n (x) be the n-fold convolution of F. 

• F € TZ~ a if F(x) = x~ a l(x) where I is a slowly varying function, i.e. 
lim^oo = 1 for all t > 0. 
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• F G ERV(a, 0) if for < a < (3 < oo and all t > 1, 

t"/ 3 < liminf £M < limsup ^1 < t~ a . 

x^oc F(X) x^oa F(x) 

• F G S if lim^oo = 2. 

Clearly = ERV(a, a). We just write F € ERV if F G ERV(a, /3) for some 
< a < /3. It can be proved that ERV C 5. The class S is called the class 
of subexponential distributions, and among others it contains the lognormal, 
loggamma as well as the heavy tailed Weibull distribution, i.e. F(x) = e~\ x l$> 
with 7 < 1. It can be shown that if F G S then i?[e tS ] = oo for all t > 0, hence 
the name subexponential distribution. 

Culminating through a series of papers [3, 33, 35, 37, 66, 67, 68], it was proved 
in [32] that when an = and Fe5, 

\ rye rt p( x \ 

4>(t,y) / — —dx, 0<t<oo. (5.1) 

r J v x 

In addition, in the same paper it was proved that if F G 7Z- a for some a > 
and JV is a renewal process independent of the Si and Wp, 

*Kt,y)~F(y) f e- ars dm(s), (5.2) 
Jo 

where m is the renewal measure of JV. In particular, if JV is the Poisson process 
then m(s) = As, hence 

V»(t >V )~ —F(y)(l - e~ art ), 0<J<co. 
or 

A generalization of (5.2) where it is only assumed that F G ERV is given in 
[15] when t = oo. Under the additional weakened assumption that claims are 
pairwise negative dependent, i.e. for i =/= j, 

P{S i <u,S j <v)<F(u)F(v), 

they show that 

/>oo 

/ F(ye rs )dm(s), 

which is easily seen to be the same as (5.2) when F G TZ- a - Another general- 
ization is provided in [78], where under some extra conditions the formula (5.1) 
is proved to hold for the absolute ruin problem when t = oo, op = and y = 
in (1.6). 

For the light tailed case with up = results are a little less explicit. Let 
K = sup{a : 7?[e a5 '] < oo}. Then it was proved in [19] and [57] that for op = 
and any e > 0, 

lim e {K - £)y ^(y) = and lim e {K+E)v ip(y) = oo. (5.3) 
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It was also shown by examples that anything can happen to lim 1 ,_ >00 e Ky ip(y). 
The result (5.3) was also proved for the absolute ruin problem with y = 0, and 
a simplified proof can be found in [78]. 

When a R > the picture is somewhat less complex. The reason for this is 
that the financial risk caused by variations in the return processs R corresponds 
to claims F £ 1Z- P with p = -=£- — 1. So when claims have lighter tails than 
this, the asymptotics is dominated by R. Various results in this direction appear 
in [20, 34, 44]. The most precise results for the model studied here are those 
in [26]. There it is assumed that ap = 0, but due to the light tailed effect of 
Wp, the result is undoubtedly valid for ap > as well. To present the results, 
remember from the beginning of Section 2 that ip(y) = 1 when p < 0, so assume 
that p > and that E[S] < oo. 

1. If for some e > 0, E[S p+£ ] < oo, then lim^oo y p ip(y) = c for some (for 
all practical purposes) unknown c. 

2. If F £ lZ- a where a < p, then 

a R a(p - a) 

3. If F G TZ- P with F(x) = x~n{x) and M < oo where M = f™ l -^-dx, 
then ip(y) is asymptotically as in Case 1 above. If M = oo then 

v 



<?rP Ji x 



In [2] similar results, but somewhat less explicit, are proved when claim intervals 
are either Erlang distributed or a mixture of two exponential distributions. 
Asymptotics when P and R are independent Levy processes can be found in 
[48]. 

There is a considerable literature devoted to the asymptotics of (3.2). In 
[43, 44] asymptotics for ip(y) and ip(nlogy,n) for n fixed is studied, mostly in 
a large deviations context. In [64, 65] asymptotics for ip(n,y) with n fixed and 
y increasing is provided. This is done for a variety of different assumptions on 
the tail behaviour of the claims, and their results are similar to those reported 
above for ip(y) in the time continuous case. A generalization of the asymptotics 
for ip(y), allowing for dependence in the {£„} process is provided in [16], where 
it is shown that ip(y) decreases at a power rate under rather general Markovian 
assumptions. 

A slightly different kind of result in the time discrete case is in [64, 65] and 
generalized in [1 1] where it is proved that when F G S, under some rather mild 
conditions 

n / k \ 

i/;(n,y)~Y l P (vI[Zi<-y . 

k=l \ i=l / 

where rj and the £j are all independent and defined as in (3.2). For more asymp- 
totic results the reader should consult the references. 
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6. Inequalities 

Several inequalities for the ruin probability exist, but the problem is that ei- 
ther the bounds are too conservative, or that they are difficult to compute. An 
example that follows directly from (2.4) is 

H{-y) < 4>(v) < 

and when F has a decreasing failure rate, the right side can be strengthened to 
■ The problem here is of course that the distribution function H is not 

known. 

An upper bound that is easier to compute, but restricted to the case with 
<yp = or = and light tailed claims is given in [5]. To explain, assume that the 
equation 

A (E\e^ s \ - l) - j(u)(p + ru) = 
has a positive solution 7*(m) for all < u < y. Then 

Further examples of inequalities can be found in [8, 34, 37, 38, 41, 75]. 



7. Stochastic interest rates 



In this section the model is slightly changed. The risk process P is still given 
by (1.3), but instead of investing in risky assets, surplus capital is invested in a 
money account. Letting r t be the rate of return of this account, capital at time 
t now becomes 

Y t =y + P t + [ r s Y s ds. 
Jo 

Assume that r follows the stochastic differential equation 

dr t = pt(r t )dt + a(r t )dW R . t , 

where \x and a satisfy the necessary conditions for existence and uniqueness, 
and where again Wr is a Brownian motion independent of Wp. Then Y is 
no longer a Markov process, but the joint process (r, Y) is a two-dimensional 
Markov process. Therefore, the probability of ruin will depend on the state of 
both these factors, and to avoid confusion with earlier notation we write ^(r, y) 
for the infinite time ruin probability. Forgetting about existence and smoothness 
properties, straightforward use of Ito's formula gives that rp(r, y) satisfies the 
following partial integro-differential equation 

^<r 2 (r)ip rr (r, y) + -a-%i>yy(r, y) + fJ>(r)i> r (r,y) + {ry + p)tp v {r, y) 

rv . . ^ 

+ X ijj(r, y - x)dF(x) - A^(r, y) = -XF(y). 
Jo 
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Here rp r (r,y) = ^$(r, y), 4> rr (r,y) = j£r^(r,2/) and so on. In addition appro- 
priate boundary conditions are needed. 

Needless to say, this is a very complicated equation, even in the pure diffusion 
case when A = 0. Assuming that A = and that r follows the simple Vasicek 
model 

dr t = (a - br t )dt + <r R dW Rit , 
(7.1) becomes the elliptic PDE 

^(T R '4i rr (r,y) + ^(Tp4> V y(r,y) + (a - br)^ r (r,y) + (ry + p)4> y (r,y) = 0, 

on (—00,00) x (0,oo). Boundary conditons are 

$(r,0) = 1, 
lim V>(r,y) = 1, 

r— * — 00 

lim i'(r 7 y) = 0, y > 0, 

r — >oo 

lim tp(r,y) = 0. 

y—*oo 

At the time of writing no results are known for this kind of problem. 
8. Minimization of ruin probabilities 

Returning to the basic model (1.1)— (1.3), we will now assume that in additon 
to investing in the risky asset the company can also invest in a risk free asset 
with return 77, where r/ < r. We denote by at the proportion of the capital 
invested in the risky asset, so that 1 — at is the proportion invested in the risk 
free asset. With this, the analogue to (1.3) is 

dRi = {{l~at)r s +a t r)dt+a t cj R dW R .t = {r f +a t {r-rf))dt+a t a R dW Rl t. (8.1) 

Clearly, if at = a, a constant, then (1.3) and (8.1) are equivalent. 

In addition we may allow proportional reinsurance. Let bt be the fraction 
of the risk retained in the company. Then the fraction 1 — bt is covered by 
the reinsurers, and we assume that the premium rate for this reinsurance is 
(1 — b f )r]p, where rj > 1. Thus rj — 1 > is an additional charge made by 
the reinsurers. This may not always be realistic, but for the problem here it is 
necessary in order to avoid trivial solutions. With this, the analogue of (1.2) is 

/ JV t 

dP t b = pil-il-btWdt + btapdWpj-btdlY. 3 * 

\i=i 

= b t dP t -(l-b t )(v-l)pdt. 

Again if bt = b, a constant, then (1.2) and (8.2) are equivalent. 

In order to be admissible, the controls at and bt must belong to certain sets 
A and B respectively. Examples for the set A can be: 
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A 


Restrictions 


(—00, 00) 


No restrictions 


(-00,1] 


Short sale allowed, but no borrowing 


[0,oo) 


No short sale, but borrowing is allowed 


[0,1] 


No short sale and no borrowing 



There are of course other possibilities as well, and with a similar set of possi- 
bilities for the set B, we see that the number of interesting cases is very high, and 
dealing with all of them is not practical. In practice the most interesting case is 
A x B = [0, 1] x [0, 1]. However, since it is simpler, but also since it is useful to 
see how much can be gained without restrictions, most results quoted here are 
for the case A X B = (—00, 00) x (—00, 00), or the case Ax B = (—00, 00) x {1} 
when there are no reinsurance options. 

As in (1.1) we can now let 

Y t aM = y + P b + [ Y^dR", 
Jo 

and ij a ' b (y) = P(T a ' b < 00) where T a ' b = inf{i : Y t a ' b < 0}. Then define 

r r (y)= min J a >\y), 

whenever the minimum exists. If ip(y) is the ruin probability without control 
options, then clearly ip a ,b (y) < ip(y). If they exist, the minimizing policies 
will be denoted by al and b* t . Sometimes it is more practical to consider total 
amount invested in the risky asset, so we write A t = a t Y^' b . Since optimal 
controls are Markov controls, we can write al = a* (Y t a b ) and so on. If there 
are no reinsurance options, i.e. B = {1}, we just write ip a (y). 

An equivalent problem is of course to maximize the survival probability 
<l> a,b {y) = 1 - i' a ' b (y)- The Hamilton- Jacobi-Bellman (HJB) equation for this 
problem is 

max {((77 + a(r - r f ))y+p{\ - (1 - 6)7 ? ))^-0 a ' fc (y) 
aeA.bea ay 

a- I (8 ' 3) 

+ ^y°W + b 2 <? 2 p)^ a > b (y) + X(E[^ b (y bS)] Hy))j = o. 

If A = this is a pure diffusion model, and existence and uniqueness of a 
smooth solution is to some extent addressed in [7]. When A > 0, existence and 
uniqueness of a solution for the case with Cp = and F having a bounded 
density has been proved through the works [25, 27, 30, 58]. 

Before discussing specific results, note first that if 77 > and € B, with 
y > yo = (77 — 1)£, using full reinsurance and having all the capital invested 
in the risk free asset, investment income is larger than premium loss due to 
reinsurance. Therefore, for y > y , a*{y) = b*{y) = and tp a *' b *(y) = 0. This 
argument is not valid when 77 = 0. 
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Unless stated otherwise, in the following the set A = (—00, 00). 

The first results pertaining directly to this problem were obtained for the 
pure diffusion case and no reinsurance, i.e. A = and B = {1}, in [7]. When 
rj = he proved that A\ = Aq > 0, so that it is optimal to invest a fixed 
amount in the risky asset. When r/ > 0, the solution is more complicated, but 
lim^oo A*(y) = 0. Also, when 77 = 0, lim^oo e K v ip a (y) = 7 for some known 
k* > and 7 > 0, while for 77 > 0, lim^oo e Ky ip a * (y) = for all k > 0. 
Thus the asymptotics differ markedly from the case without investment control 
reported in Section 5, where we saw that the ruin probability goes to zero at 
a power rate only. The reason is of course that it is optimal to invest only a 
fixed amount in the risky asset, not a fixed proportion. Further examples for the 
diffusion model, including some that allow for reinsurance, can be found in [53]. 

From now on it is assumed that A > and that ap = 0, and also that F has 
a bounded density so that (8.3) has a unique solution. 

In the light tailed case, i.e. when Ms{k) = E[e KS ] < 00 for some k > 0, 
the results do not differ very much from those in the diffusion case. Although 
A\ is not a constant, when 77 = and B = {1}, it is proved in [21, 31] that 
hm^^oo A*(y) = Aq > and also that 

lim e K * v ^ a \y) = 7 > 

y— >oo 



where n* is the positive solution of 



\{M s (k)-1)- P k= 2 



1 r 
2ct 



R 

When rj > we would expect that lim^oo A*(y) = as in the diffusion case, 
but this has not been proved. However, numerical examples given in [39] indicate 
that this conjecture holds. 

The heavy tailed case is somewhat more complicated as the results will vary 
between subclasses. When F € 1Z- a and B = {1}, it was proved in [25] that 

lim a*(y) = —-r^—, 

y— >oo <7fj 1 + a 

r'{y) - lF(y) 

for some known constant 7. Comparing with the asymptotics of Section 5, it is 
seen that controlling the investment rate only results in a better convergence rate 
when investment risk exceeds insurance risk. In that case the control reduces the 
investment risk so that it is dominated by the insurance risk. Still with B = {1}, 
the bigger class S* C S defined by 

r F(y-x)F(x) , 
lim / yy > y > dx = 2E[S] 
y^oo J Q F{y) 

has been studied in [60] when rj = and in [27] when Tf > 0. This class is rather 
big, containing lt- a for a > 1, the lognormal distribution and the heavy tailed 
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Wcibull distribution. For this reason, the asymptotics varies within subclasses 
of S* , and again whether rf = or not. For example, with F{x) = , i.e. 
heavy tailed Weibull, it is shown in [27] that 

i> a "{y) ^ j_ 
%"(y) s/y' 

for some 7 > 0. Here tpQ (y) is the minimum ruin probability when 77 = and 
ip a (y) the same when 77 > 0. 

Finally, when B = [0, 1] and F(x) > ce~ x for some c > and < e < \, it 
is proved in [59] that for 77 = there is a 7 > and a k* > so that 

lim e K ' y iP a *- b * (y) =7 

and 

lim A*{y) = A Q > 0, 

y— too 

lim = 0. 

y— >oo 

The reason we obtain an exponential rate of ip a ' b (y) is that a smaller and 
smaller fraction of the heavy losses is retained as F t a ' b increases. 

For more about the problem discussed in this section, the reader can consult 
[29, 61]. See also [5G] for a somewhat different approach. 
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